Observation of Resonant Photon Blockade at Microwave Frequencies 
using Correlation Function Measurements 
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Creating a train of single photons and monitoring its propagation and interaction is challenging 
in most physical systems, as photons generally interact very weakly with other systems. However, 
when confining microwave frequency photons in a transmission line resonator, effective photon- 
photon interactions can be mediated by qubits embedded in the resonator. Here, we observe the 

phenomenon of photon blockade through second-order correlation function measurements. The 
experiments clearly demonstrate antibunching in a continuously pumped source of single microwave 
photons measured using microwave beam splitters, linear amplifiers, and quadrature amplitude 
detectors. We also investigate resonance fluorescence and Rayleigh scattering in MoUow-triplet-like 
spectra. 



Sources of radiation differ not only by their frequency, 
but also by the statistical properties of the emitted pho- 
tons [1]. Thermal sources emit radiation that is charac- 
terized by an enhanced probability of emitting photons 
in bunches. Coherent sources, such as a laser, emit ra- 
diation with a Poisson-distributed photon number. The 
statistics of these two sources can be explained classi- 
cally. In contrast, individual atoms emit photons one 
by one well separated in time from each other, a phe- 
nomenon for which antibunching - a unique quantum 
characteristic of the field - can be observed. 

In strongly nonlinear systems, a phenomenon known 
as photon blockade [2, 3] can be used to generate a train 
of single photons that displays antibunching. Photon 
blockade is usually realized in cavity quantum electro- 
dynamics (QED) setups. Coherent radiation at the in- 
put of a cavity coupled to an anharmonic system, such 
as a single atom, is converted into a train of single pho- 
tons in the transmitted light. The nonlincarity of the 
atom-cavity system prevents more than a single excita- 
tion of the same energy entering the cavity. Only once 
the photon has left the cavity can the system be re- 
excited, realizing a single-photon turnstile device. The 
transmitted radiation has two important characteristics: 
sub-Poissonian photon statistics and photon antibunch- 
ing. On the one hand, sub-Poissonian statistics are ex- 
perimentally demonstrated when the second order cor- 
relation function fulfills the inequality .^'•^•'(t) < 1 for 
all times t. On the other hand, photon antibunching is 
demonstrated by a rise of g^"^^ (r) with r increasing from 
to larger values while g^^\Q) < g^'^\T), as discussed in 
detail in Ref. 4. 

At optical frequencies, resonant photon blockade - cav- 
ity and atom share the same resonance frequency was 
demonstrated with a single trapped atom in an optical 
cavity [5]. These measurements suffer from adverse ef- 
fects of trapping beams, micro-motion of the atom in its 
trap, and the necessity of post-selecting data for instances 



of single-atom measurements. In the solid-state, resonant 
photon blockade was demonstrated with a quantum dot 
in a photonic crystal cavity [6]. Those experiments suf- 
fered from quantum dot blinking, and limited detector 
time resolution. Our experiments are done in the mi- 
crowave regime with a single superconducting artificial 
atom resonantly coupled to a transmission line resonator, 
realizing a cavity QED setup [7] in a circuit reaching the 
strong coupling limit [8, 9]. The artificial atom at rest, 
which is here well approximated by a two-level system, 
has a strong, fixed coupling to the resonator. In ad- 
dition, our setup benefits from high-efficiency emission 
of photons in forward direction by employing an asym- 
metric quasi-one-dimensional resonator dominated by a 
single mode resonant with the artificial atom. This is 
in contrast to the atomic case for which the m\ilti-modc 
structure of the cavity is important [5]. Also, the effec- 
tive polarization of the radiation is fixed by the boundary 
conditions enforced by the superconducting metal form- 
ing the resonator, and thus does not play a role in our 
experiments. 

In this paper wc present correlation function measure- 
ments of continuous sources of single photons, coher- 
ent, and thermal radiation in the microwave frequency 
domain. In particular, wc investigate the phenomenon 
of photon blockade both in resonance fluorescence and 
second-order correlation function measurements, display- 
ing sub-Poissonian photon statistics and antibunching. 
Photon blockade in superconducting circuits has also 
been independently studied in the dispersive regime in 
Ref. 10. 

Our experimental setup is composed of two essential 
ingredients: a photon source and a quadrature ampli- 
tude detection system from which we extract the photon 
statistics. The continuous single-photon source consists 
of a single superconducting artificial atom - a transmon 
qubit [11] with transition frequency LOa - resonantly cou- 
pled to a transmission line resonator with resonance fre- 
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FIG. 1. Rayleigh scattering and resonance fluorescence of lower Jaynes-Cummings doublet, (a) Energy level 
diagram of a resonantly coupled cavity QED system driven with amplitude flu on the ground state \gO} to lower doublet |1— ) 
transition. The MoUow-type transitions arising from the dressing of the dressed states by the strong drive are also indicated 
on the side, (b) Measured resonance fluorescence spectrum including Rayleigh scattering peak (dots) at fixed drive amplitude 
of f2ij/27r = 7.9 MHz and simulated spectrum (solid line), (c) Measured resonance fluorescence spectrum vs. (indicated) 
drive amplitude Qr/2t: (dots) and analytical spectrum (solid lines). The Rayleigh peak has been omitted in these plots, 
(d) Measured MoUow side peak frequencies f^sp vs. drive amplitude Hr (dots), linear dependence ilsp = (dashed black 
lines) and calculated frequencies f^sp (solid red lines) are shown. 



quency w^/Ztt = uJa/^Tr = 6.769 GHz. In this device the 
coherent dipole coupling strength .g/27r = 73 MHz dom- 
inates over the dissipation due to photon loss from the 
cavity at rate k/2tt 4 MHz and the qubit decay at rate 
j/2tt w 0.4 MHz. When radiation impinges on the res- 
onator input at frequency {ujr — g)/2T:, only a single pho- 
ton can enter at a time, see Fig. la. Additional photons 
are prevented from entering the resonator, as transitions 
into higher excited states are blocked due to the strong 
non-linearity of the resonantly coupled qubit-resonator 
system [12-14]. In analogy with measurements in meso- 
scopic systems, where electron transport is blocked by 
the strong Coulomb interaction in a confined structure, 
this process is called photon blockade [3]. Only once the 
photon has left the cavity can the next photon enter into 
the resonator, realizing a source of single photons. 

In order to investigate the statistical properties of our 
microwave frequency radiation source, we have realized a 
scheme for measuring photon correlation functions using 
linear detectors [15-17] instead of single-photon counters 
which are still under development in this frequency do- 
main [18, 19]. In our scheme, the radiation of the source 
is passed through an on-chip 50/50 beam splitter, then 
the signal in each output of the beam splitter is ampli- 
fied using independent phase preserving linear amplifiers 
with system noise temperature T„ = 10.6 K. Finally, 



both quadrature amplitudes of each output signal are ex- 
tracted in a heterodyne measurement similar to the one 
discussed in Ref. 16. Expectation values of field ampli- 
tude, power, first- and second-order correlation functions 
can be extracted from the instantaneous values of the 
measured quadrature amplitudes. We refer to Ref. 15 for 
a detailed theoretical discussion. 

We set up our continuous single-photon source by tun- 
ing the transmon qubit transition frequency into res- 
onance with the resonator using magnetic flux [11, 20]. 
When probing the resonator transmission with a weak 
coherent tone resulting in an average resonator photon 
number (n) <C 1, we observe a characteristic vacuum 
Rabi mode splitting [8, 14] resulting from the anhar- 
monic level structure shown schematically in Fig. la. In 
many experiments of this type, only the Rayleigh scat- 
tered (elastic and coherent) part of the transmitted am- 
plitude is detected in a heterodyne measurement with a 
small effective bandwidth of ^ 50 kHz. Here however, 
we have digitally recorded the resulting fields vs. time 
in both arms of the beam splitter with a bandwidth of 
^ 50 MHz. Instantaneous power spectra of the source 
are then calculated as the product of the Fourier trans- 
form of the time-dependent signals in each arm, which 
are subsequently averaged. Here, we observe not only 
the Rayleigh scattered radiation (narrow high- amplitude 
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peak in Fig. lb) but also the incoherently scattered res- 
onance fluorescence part of the spectrum (broad low- 
amplitude triplet in Fig. lb). The resonance fluores- 
cence spectrum is characterized by three spectral lines 
(four transitions [Fig. la] , two of which are degenerate) 
forming a Mollow triplet of a resonantly driven effective 
two-level system. The two levels are realized by the joint 
ground state |f;0) and the lower energy state of the first 
doublet |1— ) = (jgl) — |e0))/-\/2 of the Jaynes-Cummings 
ladder. The dressing of these dressed states by the drive 
field has been discussed theoretically in Ref. 2 and has 
also been experimentally investigated with superconduct- 
ing circuits considering only the Rayleigh scattered part 
of the radiation [21]. 

The full spectrum is in excellent agreement with the 
numerically calculated steady-state solution of the mas- 
ter equation taking into account two qubit levels and five 
resonator levels (solid line in Fig. lb). For this calcula- 
tion, we use the device parameters quoted above and take 
into account the finite bandwidth of our detection sys- 
tem. Also, the analytically calculated fluorescence spec- 
trum of the coherently driven effective two-level system 
(solid lines in Fig. Ic) is virtually indistinguishable from 
the master equation calculation and the data. Here we 
include dephasing and do not make approximations for 
the strength of the drive [22]. To correctly capture the 
amplitude of the coherently scattered radiation in the an- 
alytical calculation the higher doublet ]2— ) needs to be 
included. 

The frequency ilgp by which the Mollow side peaks are 
offset from the central peak is observed to depend on 
the drive amplitude 51 (Fig. Ic). Lorentzian fits to the 
triplet spectrum yield fisp, showing an approximately lin- 
ear scaling with at large drive amplitudes, see Fig. Id. 
However, figp is significantly smaller than Qn at drive 
strengths smaller than the characteristic rate of dissipa- 
tion, an effect that is accurately explained by the analyt- 
ical two-level model (see red solid lines in Fig. Id) [22]. 
Similar Mollow triplet-like structures have also been ob- 
served in strongly driven superconducting flux and charge 
qubits using different detection techniques [23-25] . 

We note that, for these measurements, the uncorre- 
lated noise added by the two independent amplifiers is 
efficiently averaged out [26] , and the residual noise offset 
- a factor of 10'^ smaller than the noise introduced by a 
single amplifier - is determined by performing a reference 
measurement where the system is left in the ground state 
and then subtracted from the data [15]. 

The experiments discussed above clearly demonstrate 
the resonance fluorescence emitted from the cavity when 
it is weakly driven on the lower Rabi resonance (cj^ — (?) . 
In this limit, photon blockade is expected to be observ- 
able in measurements of the normalized second-order cor- 
relation function g(^)(r). We extract ^'^•'(t) from a mea- 
surement of the cross-correlation of the power detected 
between the two outputs of the 50/50 beam splitter [15]. 
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FIG. 2. Correlation function measurements. 

(a) Second-order correlation function measurements g^'^^r) 
(dots) for indicated drive amplitudes ^ln and master equa- 
tion calculation with and without accounting for finite 
measurement bandwidth (gray and black lines, respectively). 

(b) g^^^r) for a thermal field with mean photon number 
(nth) ~ 1.4 in the resonator, (c) (/'^'(r) for a coherent drive 
with (ric) ~ 1. 



The constant offset due to the noise added by the ampli- 
fiers is subtracted and the correlation function is normal- 
ized to unity for times t — oo. At low drive amplitudes 
[Hb/^tt = 2.5 MHz), we observe sub-Poissonian photon 
statistics characterized by ^^^•'(t) < g'-^^ (oo), which is ful- 
filled within the experimental noise. Additionally, pho- 
ton antibunching is observed since g^'^^r) rises for t away 
from T = while g^^^'iO) is at a minimum (Fig. 2a). For 
T — >■ CX3 we note that g'^'^^ approaches a constant value, 
as expected. We observe a small overshoot of g^^^r) at 
around r = tt/^Ih = 200 ns in Fig. 2a. This indicates 
a correlation between a photon emitted at time t and a 
second photon emitted with high probability at the later 
time {t + r)/r2fl = tt at which the drive has coherently 
re-excited the coupled system. 

When increasing the drive amplitude f2fl we observe 
characteristic oscillations in the measured g*^^^ (r) exactly 
at the frequency ^Ir (Fig. 2a) while we continue to clearly 
observe the signature of antibunching at r = 0. We quan- 
titatively compare the measured data to numerical cal- 
culations of ^'^•'(t) (see black lines in Fig. 2a) based on 
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a master equation calculation using the known system 
parameters. Considering the finite bandwidth ~ 20 MHz 
of the digital filter used in the quadrature data acquisi- 
tion, wc find excellent agreement between the measured 
data and the calculations, see gray lines in Fig. 2a. The 
small residual deviations of the measured g^'^^r) from 
the simulations arc due to the noise added by the am- 
plifiers. We note that each data trace was collected over 
17 hours corresponding to approximately 5.5 x 10^*^ mea- 
sured photons and 15.75 Tbyte of analyzed quadrature 
amplitude data using fast field-programmable gate ar- 
ray based electronics [16]. The presented data clearly 
demonstrates the phenomenon of photon blockade in the 
microwave domain detected using second-order correla- 
tion function measurements. 

For reference we have also measured g'-^-' (r) when pop- 
ulating the resonator with a mean thermal photon num- 
ber (nth) ~ 1-4. The quasi-thermal field distribution was 
realized by mixing a fixed frequency microwave tone with 
a large bandwidth white noise source [27]. We clearly ob- 
serve bunching g^'^^ (0) = 2 of the thermal radiation emit- 
ted from the resonator, g^"^^ (t) approaches unity on the 
time scale of the cavity decay rate k/2tt also considering 
the finite detection bandwidth (gray line in Fig. 2b). Per- 
forming a similar experiment with a coherent source de- 
rived from a strongly attenuated commercial microwave 
generator populating the resonator with (uc) « 1.0, we 
find g^^^ir) = 1 everywhere, which is in good agreement 
with the temporal statistics of a coherent source. 

We have performed correlation function measurements 
with linear quadrature amplitude detectors in the mi- 
crowave frequency domain demonstrating photon block- 
ade in a circuit QED system. We have also shown bunch- 
ing of thermal photons and probed the second-order cor- 
relation function of coherent radiation. The techniques 
and results presented in this paper have the potential to 
inspire new work controlling the flow of photons, gener- 
ating and detecting individual photons and investigating 
single-photon effects in superconducting circuits. In par- 
ticular, the observation of photon blockade will enable fu- 
ture experimental work on photon interactions in cavity 
arrays that is actively theoretically investigated [28-32]. 

The authors would like to acknowledge fruitful discus- 
sions with Barry Sanders. This work was supported by 
the European Research Council (ERC) through a Start- 
ing Grant and by ETHZ. M.P.d.S. was supported by a 
NSERC postdoctoral fellowship. A.B. was supported by 
NSERC, CIFAR, and the Alfred P. Sloan Foundation. 



[1] D. Walls and G. Milburn, Quantum Optics (Springer Ver- 

lag, Berlin, 1994). 
[2] L. Tian and H. J. Carmichael, Phys. Rev. A 46, R6801 

(1992). 



[3] A. Imamoglu, H. Schmidt, G. Woods, and M. Deutsch, 

Phys. Rev. Lett. 79, 1467 (1997). 
[4] X. T. Zou and L. Mandel, Phys. Rev. A 41, 475 (1990). 
[5] K. M. Birnbaum, A. Boca, R. Miller, A. D. Boozer, et al.. 

Nature 436, 87 (2005). 
[6] A. Faraon, I. Fushman, D. England, N. Stoltz, et al., 

Nat. Phys. 4, 859 (2008). 
[7] S. Harochc, "Fundamental systems in quantum optics," 

(Elsevier, 1992) Chap. Cavity quantum electrodynamics, 

p. 767. 

[8] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, et al.. 

Nature 431, 162 (2004). 
[9] R. Schoelkopf and S. Girvin, Nature 451, 664 (2008). 
[10] A. J. Hoffman, S. J. Srinivasan, L. Spietz, J. Aumentado, 

et al., "Dispersive photon blockade in a superconducting 

circuit," (2010), arXiv:1008.5158 [cond-mat.mes-hall]. 
[11] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, et al., 

Phys. Rev. A 76, 042319 (2007). 
[12] I. Schuster, A. Kubanek, A. Fuhrmanek, T. Puppe, et al., 

Nat. Phys. 4, 382 (2008). 
[13] M. Hoacinz, E. M. Wcig, M. Ansniann, R. C. Bialczak, 

et al., Nature 454, 310 (2008). 
[14] J. M. Fink, M. Goppl, M. Baur, R. Bianchetti, et al.. 

Nature 454, 315 (2008). 
[15] M. P. da Silva, D. Bozyigit, A. Wallraff, and A. Blais, 

Phys. Rev. A 82, 043804 (2010). 
[16] D. Bozyigit, C. Lang, L. Steffen, J. M. Fink, et al., Nat. 

Phys. 7, 154 (2011). 
[17] E. P. Menzel, F. Deppe, M. Mariantoni, M. A. 

Araque Caballero, et al, Phys. Rev. Lett. 105, 100401 

(2010). 

[18] Y. F. Chen, D. Hover, S. Scndclbach, L. Maurcr, et al., 
"Microwave photon counter based on joscphson junc- 
tions," (2010), arXiv:1011.4329 [cond-mat.supr-con]. 

[19] G. Romero, J. J. Garci'a-Ripoll, and E. Solano, Phys. 
Rev. Lett. 102, 173602 (2009). 

[20] J. Schreier, A. Houck, J. Koch, D. I. Schuster, et al., 
Phys. Rev. B 77, 180502 (2008). 

[21] L. S. Bishop, J. M. Chow, J. Koch, A. A. Houck, et al., 
Nat. Phys. 5, 105 (2009). 

[22] H. J. Carmichael, Statistical Methods in Quantum Op- 
tics 1: Master Equations and Fokker-Planck Equations 
(Springer- Verlag, 1999). 

[23] M. Baur, S. Filipp, R. Bianchetti, J. M. Fink, et al., 
Phys. Rev. Lett. 102, 243602 (2009). 

[24] M. A. SiUanpaa, J. Li, K. Cicak, F. Altomare, et al., 
Phys. Rev. Lett. 103, 193601 (2009). 

[25] O. Astafiev, A. M. Zagoskin, A. A. Abdumalikov Jr., 
Y. A. Pashkin, et al.. Science 327, 840 (2010). 

[26] G. S. Agarwal and S. Chaturvedi, Phys. Rev. A 49, R665 
(1994). 

[27] J. M. Fink, L. Steffen, P. Studer, L. S. 
Phys. Rev. Lett. 105, 163601 (2010). 

[28] S. Schmidt, D. Gerace, A. A. Houck, G. Blatter, et al., 
Phys. Rev. B 82, 100507 (2010). 

[29] M. Hartmann, F. Brando, and M. Plenio, Laser & Pho- 
ton. Rev. 2, 527 (2008). 

[30] J. Koch and K. Le Hur, Phys. Rev. A 80, 023811 (2009). 

[31] D. G. Angelakis, M. F. Santos, and S. Bose, Phys. Rev. 
A 76, 031805 (2007). 

[32] A. D. Greentree, C. Tahan, J. H. Cole, and L. C. L. 
Hollenberg, Nat. Phys. 2, 856 (2006). 



Bishop, et al.. 



